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EFFECTS  OF  ERRORS  ON  TIIE  SIDE-LOBE  LEVEL  OF 
A  LOW-SIDE-LOBE  ARRAY  ANTENNA 


INTRODUCTION 

Low-side-lobe  array  antennas  have  received  wide  interest  in  recent  years.  Theoretically,  one  can 
design  an  array  antenna  with  any  desired  side-lobe  level.  However,  in  practice  array  errors  and  other 
imperfections  limit  the  side-lobe  level.  One  question  is  then  how  low  a  side-lobe  level  one  may  achieve 
in  practice.  This  problem  is  essentially  that  of  finding  the  effect  of  array  errors  on  the  array  radiation 
pattern.  This  effect  has  been  analyzed  extensively  in  the  literature.  Ruze  [1]  considered  the  effect  on 
the  radiation  pattern  of  random  errors  in  the  exciting  currents.  Bailin  and  Ehrlich  [2]  treated  the  physi¬ 
cal  errors  which  cause  the  random  errors  in  the  exciting  currents.  Gilbert  and  Morgan  [3]  treated  the 
effect  on  gain  of  random  geometric  errors  in  the  general  two-dimensional  aperture.  Elliot  [4]  further 
treated  the  problem  of  tolerance  for  two-dimensional  scanning  antennas.  Allen  et  al.  [5]  reviewed  the 
general  problem;  they  did  extensive  study  on  all  aspects  of  this  problem. 

In  this  report,  the  relationship  between  the  side-lobe  level  and  the  array  random  errors  will  be 
examined.  In  particular,  the  limitation  on  the  side-lobe  level  as  a  function  of  the  array  errors  will  be 
presented.  Intuitively,  one  may  see  that  this  limitation  must  be  a  function  of  the  desired  side-lobe 
level,  the  number  of  elements  in  the  array,  and  the  nature  of  these  errors.  Since  these  errors  are  gen¬ 
erally  random  in  nature,  the  results  are  in  terms  of  probability  distributions.  Allen  et  al.  [5]  showed 
that  all  array  errors  that  result  from  feed,  phase  shifters,  mechanical  location,  and  the  orientation  of 
radiating  elements  can  be  characterized  by  a  phase  error  and  an  amplitude  error  for  each  element  in  the 
array.  The  results  of  this  study,  therefore,  are  in  terms  of  these  errors.  In  the  past,  these  errors  have 
generally  been  assumed  to  be  statistically  independent.  However,  there  are  cases  for  which  errors  in 
many  elements  are  not  necessarily  independent.  For  example,  the  same  phase  and  amplitude  in  a  row 
or  column  feed  network  could  feed  to  every  element  in  a  particular  row  or  column,  or  in  the  case  of 
subarray  configuration,  the  same  error  may  propagate  to  every  element  in  a  particular  subarray.  These 
errors  are  correlated  in  these  groups  of  elements.  The  effects  of  these  correlated  errors  will  be  also  dis¬ 
cussed. 

STATISTICAL  DISTRIBUTION  OF  ARRAY  PATTERN 

For  simplification,  linear  arrays  will  be  treated  first.  The  array  pattern  of  a  linear  array  can  be 
represented  by 

P(0)  =  ^A„exp[~j(2nnd/\)(sin0  -  sin0o)l.  0) 

n 

where  0  is  the  angle  of  incidence  of  a  plane  wave  on  the  array  and  0o  is  the  beam  pointing  angle.  Ele¬ 
ment  spacing  d  is  assumed  to  be  uniform.  We  further  define 

fi  -  ( 2irt//\)(sin  0  -  sin  0O).  (2) 

Equation  (1)  then  becomes 

expfjnn).  (3) 

n 

For  radiation  in  real  space,  fi  is  constrained  so  that 

_  Ul  ^  2ir.  (4) 
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Due  to  mechanical  and  electrical  errors  in  the  array,  the  array  pattern  becomes 

£/4,,(l  +8„)  exp O'dO  exp (Jnn),  (5) 

If 

where  is  the  amplitude  error  in  percent  and  <fin  is  the  phase  error.  These  errors  vary  from  element 
to  element  and  are  random  in  nature.  For  simplification,  we  assume  that  these  errors  have  a  known 
probability  density  function.  Because  of  randomness  of  these  errors,  G(fi)  is  a  random  complex  func¬ 
tion  which  is  the  sum  of  many  random  variables.  Each  of  these  random  variables  can  be  represented  as 


-  (l+5„)  exp(j<t>„)  exp(/n*i) 


(6) 


-  Xn  +JYn. 


These  random  variables  are  independent  and  have  the  same  probability  density  function.  The 
array  pattern  is  hence  a  random  function  of  the  sum  of  many  random  variables,  so  that 

Gip.)  =  g|(/t)  +  jg2(p.) 

(7) 

n  n 


According  to  Lindenberg  and  Levy’s  central  limit  theorem  [6),  g,  (p.)  and  g2(p.)  are  asymptotically 
normal  The  means  and  variances  of  £i(/x)  and  g2(p.)  are  the  weighted  sums  of  means  and  variances  of 
Xn  and  Y„.  The  means  of  g,(/i)  and  g2(jx)  are  then  respectively 


g2(pi)  =  <J>(l)^/ln  COS 

n 

(8a) 

and 

g2(fl)  =  <J>(l)£/ln  sin  nP- 

n 

(8b) 

where  d>  \k )  is 

the  characteristic  function  of  random  variable  x,  defined  as 

=  J  g(x)  exp  (jkx)  dx, 

where  g(x)  is  the  probability  density  function  of  random  variable 
above  expression  we  have  also  assumed  that  the  amplitude  error  8„ 
is  symmetrically  illuminated,  such  that  An  =  A_„,  then  g 2(ju)  =0. 
are,  respectively, 

x.  Furthermore,  in  deriving  the 
has  zero  mean.  If  the  linear  array 
The  variances  of  g[(jx)  and  g2(jx) 

o-j2  =  l/2£  Aj(A  +  B  cos2n/i) 

n 

(9a) 

and 

crl'**  1/2  A} (A  —  B  coslrifi). 

If 

(9b) 

where 

A  -  1  +tr,2-*1(l) 

(9c) 

and 

B~(  1  +  o-g2)<t>(2)  -  d>2(l). 

(9d) 

The  quantity  <r82  is  the  variance  of  the  amplitude  error  8.  The  covariance  of  giip,)  and  g2(p.)  is 

<r,2  -  l/2£  A*B  sin  2 ngi. 


(10) 
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For  the  symmetrically  illuminated  array,  cr12  is  zero.  Derivations  of  these  variances  are  included  in 
Appendix  A.  The  joint  density  function  of  the  complex  variable  Pip.)  is 


P[g\ip).  g2ip )1  = 


2lTyJ(T\(T2 


exp 


(gi~ £i>2  _  ui 
2<t\  2<t  2 


(11) 


The  probability  density  function  P{gi,g2)  is  a  generalized  noncentral  chi-square  distribution  with 
two  degrees  of  freedom.  When  g,  —  0  and  cr,  “  cr2,  this  becomes  a  Rayleigh  distribution.  For  the 
case  ^  0  and  cr ,  —  cr  2,  it  is  sometimes  referred  to  as  a  Rician  distribution.  Since  the  chi-square  dis¬ 
tribution  cannot  be  evaluated  in  a  straightforward  way,  we  shall  attempt  to  determine  if  it  can  be 
approximated. 


The  variances  o'!  and  cr2,  as  shown  in  Eqs.  (9a)  and  (9b),  consist  of  two  parts.  The  first  part  is 


P\  =  1/2  A  £  Aj. 

n 


(12a) 


and  the  second  part  is 

p2  =  1/2  B  X  'in  cos  2np.  (12b) 

n 


The  first  part  is  not  a  function  of  p,  but  the  second  part  is.  Since  p2  is  the  sum  of  cosine  functions 
except  at  regions  in  the  vicinity  of  p  -  kir ,  p,  is  much  greater  than  p2-  An  example  is  shown  in  Fig.  1, 
on  whicn  two  sets  of  curves  are  shown.  Curve  1  shows  the  cr  i  and  cr2  values  as  a  function  of  p  for  the 
case  of  a  20-element,  30-dB  Chebyshtv  array.  The  amplitude  error  has  a  variance  of  10%,  and  the 
phase  error  has  a  normal  distribution  with  zero  mean  and  a  5°  variance.  One  may  see  that  the  vari¬ 
ances  are  independent  of p  and  cr,  -  cr2  except  in  the  regions  at p  ~  it  and/u.  -  27r.  Curve  2  on  this 
figure  shows  the  case  of  an  80-element  Chebyshev  array  with  50-dB  side  lobe  design,  an  amplitude 
error  of  0.005,  and  a  phase  error  of  2°  rms.  Note  that  cr ,  and  cr  2  of  curve  2  are  considerably  smaller 
than  those  of  curve  1.  This  is  because  there  are  more  elements  and  a  smaller  error  in  the  case  of  curve 
2.  This  will  be  discussed  in  more  detail  later.  From  this,  one  may  conclude  that  in  the  side-lobe  region 
where  p  &  tt  one  may  assume  that  cr ,  —  <r2  —  px.  The  chi-square  probability  density  of  Eq.  (11)  then 
becomes  Rician,  and 


U(n) 

Fi*.  1  —  Values  of  <r,  and  cr2  (case  1:  20  element  Chebyshev,  30  dB, 
<Tj  —  0.1,  <r+  —  5°;  case  2:  80  element  Chebyshev,  50  dB,  < rt  —  0.00$, 
<r*-2T 


3 
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P(R)“  -4  exp  [-  (R  -  *,)2]  /0  Py-  ,  (13) 

where  R  -  y/g\  +  g2  is  the  amplitude  of  the  radiation  pattern  and  10  (z)  is  the  modified  Bessel  func¬ 
tion  of  zero  order. 

Radiation  in  regions  in  the  vicinity  of  n  *=  kir,  which  contain  the  main-beam  and  grating-lobe 
regions  and  a  small  portion  of  the  side-lobe  region,  has  essentially  a  chi-square  statistical  distribution  as 
shown  in  Eq.  (11). 

Before  further  discussing  the  radiation  level  statistic,  let  us  examine  a  \  and  cr 2  more  carefully. 
VALUES  OF  cr,  AND  a2 


In  the  vicinity  of  n  =  ki r,  <t\  and  <r2  can  be  approximated  by 


trf  =  1/2  (A  +  B)  £  A* 

n 

(14a) 

and 

err  ^  1/2  (A  -  *)  £  A2. 

fj 

(14b) 

and  in  other  sidelobe  regions  one  has 

<rf  =  cr?  *=  A  £  A2. 

(14c) 

Each  <r  consists  of  two  parts.  The  first  part  is  £  A „2,  which  is  a  function  of  the  number  of  elements  in 

ft 

the  array  and  the  illumination  function  of  the  array.  The  second  part  is  A,  A  -  B,  or  A  +  B  and  is  a 
function  of  the  parameters  A  and  B  given  in  Eqs.  (9c)  and  (9 d),  which  are  determined  by  random  array 
errors. 

In  finding  the  first  part  nf  cr ,  for  convenience  of  comparison,  let  us  normalize  the  array  illumina¬ 
tion  function  A„  in  such  a  w,  rat 

^4,-1.  (15) 

This  implies  that  at  the  peak  of  the  main  beam  the  radiated  field  has  unit  strength  (or  zero  dB).  The 
summation  of  A 2  is  then  always  less  than  unity.  In  the  case  of  a  uniformly  illuminated  array, 

(16) 

n 

where  N  is  the  total  number  of  antenna  elements.  In  Appendix  B  it  is  shown  that  no  matter  how  the 
illumination  function  is  formed  the  1/ N  factor  is  its  lower  bound;  in  other  words, 

I  ^  (17) 

n 

Figure  2  shows  some  examples  of  £  A*  for  a  Chebyshev  array.  It  is  evident  that  in  all  the 

n 

different  designs  the  £  A}  is  close  to  the  factor  1/  N  and  it  is  always  bounded  by  this  factor.  This  is 

n 

useful  for  the  estimation  of  <r .  One  may  replace  X-4,2  by  the  inverse  of  the  number  of  elements  in  the 
array  for  a  first-order  estimation. 


-w 
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Fig.  2  —  Values  of  £  f°r  »  Chebyshev  array  with  different 

n 

side-Jobe  levels 


In  finding  the  second  part  of  cr ,  we  consider  the.  phase  error.  The  probability  density  function  of 
the  phase  error  <f>  is  an  even  function  with  zero  mean;  hence  4>(Ar)  is  real.  Furthermore,  in  most  cases 
of  interest,  the  probability  function  is  most  likely  concentrated  in  a  very  small  angle  region  (less  than 
±  n/2).  In  this  case  d>(l)  and  <J>(2)  would  be  positive  and 

<*>(*)  <  •( 0)-  1. 

When  the  phase  error  probability  density  is  an  impulse  function  centered  at  zero  degrees,  both 
<t»  ( 1 )  and  (2)  are  equal  to  unity.  Under  this  condition  A  =*  B  =  tr62,  the  variance  of  the  amplitude 
error.  As  the  phase  error  increases,  the  probability  density  function  p(<i>)  spreads  out  and  <t>(k)  is  less 
than  unity;  parameter  A  then  increases  monotonicalfy  as  the  phase  error  increases.  This  is  shown  in 
Figs.  3a  and  3b,  Figure  3a  shows  the  case  for  which  the  phase  error  density  function  is  normal;  Fig.  3b 
shows  the  case  for  a  uniform  density  function.  In  the  side-lobe  region,  o  =  A^A„\  one  may  therefore 

n 

estimate  the  value  very  easily  by  multiplying  this  factor  by  1/A.  Some  typical  values  of  the  parameter  B 
are  shown  in  Figs.  4a  and  4b  for  normal  and  uniform  angle  error  density  functions,  respectively. 

STATISTICAL  DISTRIBUTION  OF  SIDE-LOBE  LEVEL 

In  the  main-beam  region  where  n  ~  0  and  o-j  ^  oj,  the  radiation  level  has  a  noncentral  chi- 
square  density  function.  For  this  type  of  density  function  most  of  the  probability  mass  concentrates 
within  an  ellipse  with  major  and  minor  semiaxes  of  two  to  three  times  cr|  and  o2,  centered  at  the  mean 
value  g\.  Furthermore, 

g )  »  cr i  or  rrj- 

Therefore,  for  practical  purposes  one  may  assume  that  the  radiation  amplitude  (/?o  =  V W  +  g\ ),  is 
equal  to  g,  with  a  probability  of  unity, 

Rq~  g]  =  4>(1)  X  ^n- 

n 

Since  the  illumination  function  A„  is  normalized,  one  finds 

Rq  =  $>(1).  (18) 
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Fig.  4(b)  —  Values  of  B:  angle  error  has  a  uniform  distribution 


In  the  side-lobe  region,  there  are  cases  in  which  ^  -  ir  and  in  which  fi  ^  n.  Since  the  region  in 
the  vicinity  of  =  7r  is  limited  to  about  one  beamwidth,  we  shall  first  examine  the  case  for  fj,  ^  n. 

The  probability  density  function  in  the  side-lobe  region  where  fi  ^  tr  has  a  Rician  distribution. 
The  side-lobe  level  is  defined  as  the  ratio  of  the  main-beam  level  to  the  side-lobe  level.  Its  probability 
function  is  the  joint  probability  function  of  R0  and  R.  Since  we  showed  earlier  that  R0  =  4>(1)  with  a 
probability  of  unity,  the  probability  density  function  of  the  ratio  of  R  to  R0  is  the  same  as  that  of  R 
with  a  scale  factored);  that  is,  R'  =  /?/<J>(l),  and  the  density  function  of  R'  is  then 


7 
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where 


P(R')~  exp  ({(/TOO)]2-  ^!2)/2cr2) /0 


(19) 


and 


<t  =  cr  j  =  a  2 


Hi  «  <t>(l)  cos  "P 

n 


For  convenience,  we  shall  normalize  this  R '  in  such  a  way  that 

5  | 

<T*\ 

where 

s  =  R'/^An  cos  nfj. 


P(S)  - 


■  exp 


S2  +  1 


2cr’ 


(20) 


(21a) 


and 

cr’  =  <r  ;| ]T  A„  COS  .  (21b) 

We  notice  that  £  .4„  cos  ///a  represents  the  side-lobe  level  at  the  angle  /x  when  there  is  no  error 

n 

present.  Therefore,  both  S  and  a'  in  this  equation  are  measured  in  terms  of  the  designed  side-lobe 
level.  This  is  more  convenient  to  use  than  are  values  in  terms  of  <r  or 


The  cumulative  probability  of  5  being  less  than  S,  is  then 

p(s<sL)~  CL  -4 

L  *'0  „ 


exp 


S 2  +  1 

/ 

s 

2cr  '2 

cr'2 

dS. 


(22) 


A  family  of  such  curves  with  cr'  as  parameter  is  shown  in  Fig.  5.  Each  of  these  curves  presents 
the  cumulative  probability  that  S  (in  terms  of  designed  side  lobe)  is  less  than  or  equal  to  a  level  SL  for 
a  given  cr '.  Since  this  curve  is  presented  in  such  a  way  that  it  is  not  a  function  of  the  angle  n ,  these 
curves  apply  to  all  points  in  the  side-lobe  region.  Secondly,  these  curves  are  normalized  with  respect  to 
the  ideal  side-lobe  level.  It  repr'  sents  the  probability  of  the  deviation  of  side-lobe  level  from  the 
designed  value.  Although  they  are  not  presented  explicity  as  a  function  of  m  ,  they  are  related  to  the 
side-lobe  level.  For  example,  at  the  peak  of  a  side  lobe,  the  normalized  cr’  may  be  only  equal  to  0.1, 
but  at  a  point  where  the  side-lobe  level  may  be  10  times  smaller,  the  normalized  < r  ’  then  becomes  10 
times  larger.  One  can  see  the  difference  in  the  probability  distribution  for  these  two  cases.  This  set  of 
curves  is  universal.  It  applies  to  arrays  with  different  illumination  designs,  different  sizes,  and  different 
errors. 

With  the  aid  of  this  plot,  one  may  easily  determine  the  required  error  tolerance  to  achieve  a 
desired  side-lobe  level.  For  example,  one  may  wish  to  design  an  array  having  a  50-dB  side-lobe,  with 
probability  of  90%  that  the  side-lobe  level  will  not  exceed  the  designed  level  by  more  than  30"  The 
curve  for  <t'  =  0.2  satisfies  this  condition,  because  for  S  =  1.3  the  cumulative  probability  is  90%. 
Since  cr’  =  0.2  is  equivalent  to-14dB.  the  required  a:  is  approximately  64  dB  <50+14,  Suppose  that 
the  array  has  100  elements.  This  value  of  A' gives  at  best  20  dB;  one  therefore  needs  an  error  level  that 
would  yield  no  less  than  44  dB.  From  Fig.  (3a),  one  sees  that  an  amplitude  error  of  0.5%  (1  cr)  and  an 
angle  error  of  1.5°  (I  cr)  will  satisfy  this  requirement. 
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Fig.  5  —  Cumulative  probability  distribution  of  normalized  side-lobe  level 
(normalized  to  the  designed  value) 


The  validity  of  the  curves  shown  in  Fig.  5  has  been  checked  by  means  of  computer  simulation.  In 
this  simulation,  we  first  computed  the  probability  curve  as  shown  in  Fig.  S  for  an  array  of  52  elements 
with  a  50-dB  side-lobe  design  with  given  phase  and  amplitude  statistical  error  distributions.  The  results 
are  plotted  in  Fig.  6.  Next,  we  computed  the  side  lobe  level  of  this  array  about  10,000  times.  Each 
time  the  actual  phase  and  amplitude  for  each  element  of  the  array  was  generated  by  the  random - 
number  generator  according  to  the  prescribed  distribution.  These  phase  and  amplitude  errors  were  then 
added  to  each  element  in  the  array.  Pattern  values  were  computed  in  the  side-lobe  region  and  normal¬ 
ized  to  the  designed  side-lobe  level.  The  statistical  distribution  of  these  computations  is  plotted  in  Fig. 
6  with  crosses,  and  the  theoretical  curve  is  also  plotted.  One  may  see  that  these  two  results  closely 
coincide. 

SIDE-LOBE  REGION  FOR  /*  =  it 

In  the  side-lobe  region  where  n  =  ir,  cr\  <r2,  and  the  probability  density  function  of  the  radia¬ 
tion  pattern  becomes  a  noncentral  chi-square  distribution  similar  to  that  in  the  main-beam  region.  The 
relationship  that  >>  a\  and  cr2  does  not  hold;  therefore,  the  approximation  used  in  the  main-beam 
region  cannot  be  applied.  Because  the  region  in  the  vicinity  of  ^  =  n  is  very  small  (about  one 
beamwidth),  it  has  been  ignored  in  the  past.  However,  it  is  worthwhile  to  investigate  the  probability 
that  an  undesired  high  side-lobe  may  appear  in  this  area. 

It  was  pointed  out  earlier  that  most  of  the  probability  mass  of  a  noncentral  chi-square  distribution 
concentrates  within  an  ellipse  with  major  and  minor  semiaxes  of  two  times  <7,  and  cr2,  centered  at  the 
mean  value  gt.  This  situation  is  depicted  in  Fig.  7.  When  <7,  =  <t2,  the  maximum  value  of  R  is 

(23) 

and  when  a-  2  <  <  cr  i 

Si  +  2or-  (24) 

This  is  the  maximum  side-lobe  level,  with  a  very  high  probability  that  no  side-lobe  level  will 
exceed  this  value.  When  g[  >  a-^  or  o-2,  the  difference  between  values  of  Rmtx  in  Eqs.  (23)  and  (24) 
is  very  small.  Thus  for  cases  which  have  a  very  high  probability  of  a  small  deviation  of  the  side  lobe 
from  the  designed  value,  the  approximation  of  the  noncritical  chi-square  distribution  with  a  Rician  dis¬ 
tribution  is  acceptable. 
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In  this  region,  from  Eqs.  (9a)  and  (9b), 

W*  \u+B)'2,A*  (25a) 

2 

and 

<T1~  \(A-B)Y,An-  (25b) 

n 

where  A  and  B  are  functions  of  the  characteristic  function  of  phase  error  and  the  standard  deviation  of 
amplitude  error  as  shown  in  Eqs.  (9c)  and  (9d).  For  a  very  small  phase  error  both  <I>(1)  and  d>(2)  are 
close  to  unity,  hence  A  =  B,  cr\  =  2A,  and  cr2  =  0.  As  the  phase  error  increases,  A  increases  and  B 
decreases.  One  may  see  this  from  Figs.  3a  and  4a  (or  Figs.  3b  and  4b).  In  this  case,  one  may  assume 
that  B  =  0  and  <r  \  —  <t2  =  A.  The  maximum  side-lobe  levels  with  a  high  probability  for  these  two 
cases  are,  respectively, 

7?max  =  g\  +  4A  (26a) 

and 

*  max  =  g\  +  2 A.  (26b) 


One  may  see  this  intuitively.  When  the  phase  error  is  zero,  amplitude  errors  have  a  good  chance 
to  line  up  and  induce  a  higher  side-lobe  level.  On  the  other  hand,  if  phase  errors  are  introduced  the 
chance  of  all  amplitude  errors  to  be  lined  up  in  the  same  direction  is  greatly  reduced,  hence  the  R^x 
value  has  a  higher  probability  of  being  smaller.  If  the  phase  error  is  further  increased,  eventually 
A  ~  -B.  In  this  case  a-\  =  0  and  ar2  =  2 A,  and  under  this  condition 

K™,  =  V«f  +  (4  A)2.  (27) 

This  f?max  value  is  slightly  greater  than  the  one  in  Eq.  (26b).  This  is  the  case  in  which  phase 
errors  are  so  great  that  the  imaging  component  dominates  and  it  has  a  good  chance  to  be  lined  up  and 
introduces  a  higher  amplitude  error.  However,  if  gj  »  A ,  the  error  introduced  by  the  approximation 
o- 1  =  cr2  —  A  is  small.  This  is  shown  in  Fig.  8.  Curve  1  is  the  cumulative  probability  of  the  normalized 
amplitude  for  cr\  =  <r2  with  the  normalized  cr,  =  0.1;  curve  2  is  for  the  case  <7,  =  0.199  and 
<r2  =  0.0001;  and  curve  3  is  for  the  case  o-,  =  0.0001  and  c r2  =  0.199.  The  difference  in  the  probability 
distribution  of  these  three  cases  is  negligible  for  practical  purposes. 

PLANAR  ARRAYS  WITH  CORRELATED  ERRORS 

For  a  planar  array,  if  errors  in  each  element  are  independent,  the  results  of  the  linear  array  can  be 
applied  directly.  The  array  pattern  of  a  planar  array  with  independent  errors  can  be  represented  by 


Gifi.v)  =  II  ^(1  +  exp  (J<t>nm)  exp  [jimii  +  nv)],  (28) 

n  m 

where 

H  =  (2ndx/\)  (sin  9  cos  <b  —  sin  0q  cos  <£o)- 
v  =  (7ndJ\)  (sin  9  sin  <f>  -  sin  90  sin  </> 0). 

and  dnm  and  d>nm  are  the  amplitude  and  phase  errors  respectively.  It  can  be  shown  that  the  a\  and  cr2 
of  this  case  are  given  by 
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J 

3.0 

Fig.  8  —  Probability  distribution  at  region  fj.  =  rr,  er  =0.1 


W-  Vi  II  l Cm  +  Dm 

n  m 

cos  (2 mu  +  2 nv)]. 

(29a) 

0-2  ”  Vi  II  [Cnm  -  D„m 

n  m 

cos  (2 m/i  +  2 nv)], 

(29b) 

where 

Cnm-  Il+ffi,-**,  0)]  Aim 

(29c) 

and 

Dnm  -  10  +<ri.)  (2)  -  (DMi. 

(29d) 

and 

®T2  -  Vi  II  1(1  +  cr„2m)  <t>„m(2) 

-  4>n2m  sin  (2 mu  +  2 nv). 

(30) 

n  m 


Comparing  these  equations  with  Eqs.  (9a),  (9b),  and  (10)  one  may  see  that  they  are  almost  identical. 
Therefore,  all  results  obtained  for  the  linear  array  can  be  applied  directly.  Since  a  planar  array  usually 
has  more  elements  than  a  linear  array,  <t\  and  cr2  in  general  are  smaller.  It  therefore  can  tolerate  larger 
errors  with  less  degradation.  Unfortunately,  in  most  cases  errors  in  each  element  are  not  independent. 
They  may  be  correlated.  For  example,  a  planar  array  may  be  fed  by  rows  and  columns.  The  same  error 
from  the  feed  network  may  appear  in  every  element  in  an  entire  row  (or  column)  or  in  many  cases  all 
elements  of  a  subarray  may  contain  the  same  error  due  to  mechanical  reasons.  To  examine  this  effect, 
we  shall  first  examine  a  case  for  which  the  same  error  appears  in  every  element  of  a  row  (or  column)  in 
an  array.  Besides  this  correlated  error,  there  are  also  independent  errors  in  every  element.  The  array 
pattern  can  then  be  assumed  to  have  the  following  form: 

Gifi.v)  =  £0  4-  8„)  exp  m  (1  +  8„m)  exp  exp  \j(mn  +  ni>)).  (31) 

n  m 

where  5*  and  <t>„  are,  respectively,  the  amplitude  and  phase  errors  which  appear  in  the  nth  row.  Appen¬ 
dix  C  shows  that 

<t}  -  E  +  h  (32a) 


12 


NRL  REPORT  8485 


and 

o-22  =  E  -  F. 


where 

E  =  VEL  Kl  +  (1  +  "D  ~  <Pn2  0)  *L  (1)1  AL 

n  m 

+  (1  +  <T2)  -  <t>2  (1)]  <t>*„  (1)  YLAnmAns  COS  (ffl  -  s)fX 

n  ms 

and 

F-VtZZ  f(l  +  O  <»  +  (2)4>nm  (2)  -  <nfl2(l)4>2m(l)] 

n  m 

■  A2„  cos  (2m/z  +  2n«/)  +  >/4  ^  +  ^  *»^2)  -  <J>2(1)] 

n  m  s 

■  4>2m  (1)  cos  [(m  +  s)  /LI  +  2»k], 


(32b) 


(32c) 


(3  2d) 


In  the  derivation  of  Eqs.  (32),  it  is  assumed  that  the  amplitude  errors  8„  and  Snm  have  zero  mean 
and  their  respective  variances  are  cr2  and  cr  2m.  In  the  formulation,  we  also  assumed  that  the  amplitude 
error  consists  of  two  levels  8„  and  8„m,  and  the  total  error  is  (1  +  8„)  (1  +  8„m).  The  mean  of  this 
error  is  zero;  however,  the  composite  variance  is  given  by 

(1  +  a2)  (1  +  a2m)  =  1  +  ar 2  +  <r2m  +  cr2a-2m.  (33) 

Let  the  total  amplitude  variance  a  be  given  by 

<r2  “  or n  +  or 2m  +  Or2  cr2m.  (34) 

If  both  <rn  and  <rnm  are  small,  the  total  variance  a-2  can  be  viewed  as  the  sum  of  the  individual  vari¬ 
ances  <r2  and  a2m.  Let  the  total  phase  error  4>  be  the  sum  of  the  two  phase  errors  <t>„  and  4>nm,  then 
the  characteristic  function  of  <J>  is 

(1)^(1).  (35) 

Equations  (32a)  and  (32b)  consist  of  two  parts,  similar  to  Eqs.  (9a)  and  (9b)  in  the  case  of  a 
linear  array.  The  terms  which  involve  cos  (2 mu  +  2 nv)  and  cos  t(m  +  s)  n  +  2nv]  can  be  ignored 
except  perhaps  in  regions  very  close  to  the  main  beam  and  where  2 mfi  +  2 nv  and  (m  +  sV  +  2 nv  are 
integer  multiples  of  2rr.  These  regions  are  only  one  beamwidth  in  extent,  and  their  effect  may  be 
ignored  for  the  same  reason  that  was  used  in  the  case  of  a  linear  array.  We  hence  have 

o- 12  -  o- 22  -  1 h  ££  [1  +  <r2  -  4>2(l)]i4£,  +  'h  £  [(1  +  cr  2  -  <I>2(1)]  | G„  (/li)I2.  (36) 

n  m  n 

where 

I G„ |2-  <t>2„  (l)L5Xm/l«  cos  ( m  -  s)n,  (37) 

m  s 

and  cr„ ,  c rnm,  <!>„,  and  <t>nm  are  replaced  by  a-  and  <t>  as  shown  in  Eqs.  (34)  and  (35). 

The  mean  radiation  pattern  of  the  linear  array  of  the  wth  row  is 

G(p.)  =  4>nm(l)Z'4„m  exp  (Jmn),  (38) 

m 

therefore  |G(*i)|2  is  its  power  pattern.  The  <ri  and  o-2  consist  of  two  parts.  The  first  part,  1/2  ££ 

n  m 

[1  4-  o-2  -  <I>2(1)]  A2m,  is  identical  to  that  of  an  array  with  uncorrelated  errors  when  both  the  correlated 
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and  uncorrelated  errors  in  each  element  are  taken  into  account.  The  second  portion,  '/?  £ 

ft 

[(1  +  <r2  -  <t>2(D]  |C7(ju)!J,  has  a  value  similar  to  that  of  a  linear  array  with  the  illumination  weights 
A replaced  by  the  power  radiation  pattern  function  |G(/t)|2.  It  is  interesting  to  note  that  whence  =  0, 
|G(^i)P  has  its  maximum  value.  In  other  words,  at  /u.  =  0  the  correlated  errors  have  the  strongest 
effect  on  cr  and  the  array  patterns  would  most  probably  deteriorate  more  at  /x  —  0  in  the  n .  v  domain. 
For  example,  a  column  and  row  fed  array,  as  discussed  in  this  section,  has  rows  in  the  x  direction  and 
columns  in  the  y  direction.  Elements  in  each  column  have  a  common  error.  If  the  array  beam  is 
steered  at  the  broadside  ( 60  =  0),  then  n  =  sin  0  cos  <£,  and  the  maximum  degradation  will  occur  in  a 
plane  along  the  y  direction  (the  E  plane),  when  <f>  -  90°  However,  if  elements  in  each  row  have  a 
common  error,  the  worst  degradation  may  occur  at  v  =  0,  along  the  x  direction  (the  H  plane).  Under 
this  condition. 


|Gm(^)l2  =  I(I^m|2. 

(39) 

n  \  m  ) 

If  the  illumination  coefficient  Anm  values  are  normalized  such  that 

1. 

n  m 

(40) 

one  can  show  that 

(41) 

where  N  is  the  total  number  of  rows  which  have  correlated  errors.  One  may  hence  conclude  that: 


•  cri  =  <r2. 

•  Both  cr  ,  and  <j  2  consist  of  two  parts,  with  the  first  part  due  to  total  errors  (including  both 
correlated  and  uncorrelated)  at  each  array  element  (cr,, )  and  the  second  part  due  to  the  correlated  error 
fed  to  each  row  (or  column)  (<rc).  The  <r  is  the  root  sum  square  <RSS)  of  these  two  parts. 

•  Both  and  <r  c  are  functions  of  both  the  array  element  errors  and  the  illuminations.  The 
effects  of  amplitude  and  phase  errors  on  the  cr  value  are  identical  to  that  of  the  case  of  a  linear  array. 
Curves  shown  on  Figs.  3a  and  3b  can  be  used  to  estimate  this  value.  The  array  illumination  can  be 
approximated  by  the  number  of  elements  in  the  array  (for  cr„)  and  the  number  of  rows  (foro-c)  for  the 
correlated  errors. 

The  covariance  o-i2  between  real  and  imaginary  parts  of  the  array  pattern  is 
o-|2  =  'h  ££  [(1  +  cr2)  (1  +  <r (2)<t>„m (2)  -  ^(Od^a)]^  sin  (2 mu  +  2 nv) 

n  m 

+  *  III  [(1  +  o-2)d>„( 2)  -  <t>}(l)]<&2m'4„mAms  sin  Km  +  s)fi  +  2nv].  (42) 

n  m  s 

When  the  array  is  symmetrically  illuminated  and  the  phase  center  of  the  array  is  taken  at  the  center  of 
the  array,  one  may  show  that  cr  12  -  0.  Therefore,  the  distribution  of  array  patterns  of  a  planar  array  is 
Rician  as  shown  in  Fig.  5.  The  value  of  cr  is  the  RSS  of  crM  and  <r  c.  The  above  results  can  be  applied 
to  an  array  which  is  fed  by  subarrays  and  the  correlated  error  appears  in  each  subarray.  In  this  case, 
one  has  to  replace  the  parameters  fi  and  v  with  functions  of  m  and  v . 
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CONCLUSIONS 

In  this  report,  we  have  shown  the  following: 

1 .  The  amplitude  distribution  of  the  radiation  pattern  for  both  a  linear  array  and  a  planar  array  is 
Rician.  A  set  of  universal  curves  for  such  distributions  is  shown  on  Fig.  5.  The  curves  are  presented 
for  different  o-  values  and  are  normalized  to  the  mean  pattern  value.  They  thus  read  directly  in  terms 
of  cumulative  probability  of  the  degradation  of  the  side-lobe  level. 

2.  The  cr  value  for  both  linear  arrays  and  planar  arrays  has  a  similar  form,  which  is  a  function  of 
errors  and  array  illuminations.  The  effect  of  errors  is  plotted  in  Fig.  3  for  different  error  distributions. 
Figure  3  also  shows  that  the  illumination  function  of  the  array  can  be  approximated  by  the  number  of 
elements. 

3.  For  a  planar  array  with  correlated  errors  the  variance  cr  is  the  RSS  of  crM  and  a c.  The  variance 
<rM  takes  into  account  the  error  in  each  element  (including  both  correlated  and  uncorrelated  errors)  and 
arc  is  due  to  the  contribution  of  correlated  error  in  a  subarray.  Both  crM  and  vc  are  functions  of  the 
errors  and  the  illuminations.  However,  <yf  is  also  a  function  of  the  subarray  pattern  in  which  the  errors 
are  correlated. 
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Appendix  A 

DERIVATION  OF  VARIANCES  FOR  A  LINEAR  ARRAY 


Define  a  complex  random  variable 


z  -  x  +jy\ 

(Al) 

then 

E[\z  -  £(z)|2)  =  <rj  +  <rj. 

(A2) 

where 

*}-  £{[*-  E(x)]2). 

(A3) 

o-l  -=  £{Lv  -  fCv)]2}. 

(A4) 

and 

E[[z  —  £(z)]2}  =  W2  —  (r|  +  2j(T\i, 

(A5) 

where 

(t n  =  £{U  -  £(x)]  [y  -  £(y)l). 

(A6) 

Solving  Eqs.  (A5)  and  (A2),  one  gets 

2or 2  =  £{|z  -  £(z)!2}  +  Re  £{[z  -  £(z)]2}, 

(A7) 

2<t\  =  £{|z  —  £(z)|2}  -  Re  £{ [z  -  £(z)]2l, 

(A8) 

and 

2cr  ,2  -  Im  £  {{z  -  £(z)]2}. 

(A9) 

For  a  linear  antenna,  the  pattern  function  is 

GW)  -  £ A„  (1  +  5„)  eJ *•  e>^, 

n 

(A10) 

£(G)  =  <P(1)  ZA „  e^, 

(All) 

|£(G) |2  =  <J>2(1)  JX-Mm  e>Kn-m)», 

n  m 

(A12) 

[£(G)]2  =  <1>2(1)  + 

n  m 

(A13) 

£{|Gl2}-  £  A„Am  <I>2(1)  e>(n~m)»  +  £/t„2  (1  +o-2), 

(A  14) 

n^m 

where  <r2 

is  the  variance  of  the  amplitude  error  8„, 

£{|G  -  £(G)|2|  =  £{|z|2  -  |£(z)i2) 

n  n 

(A15) 

£{(G)2)  -  £  *2d)  eJ<n  +  m)>l  4-  ZAn  (1  +  a2)  <t>(r)  ej2n/i. 

n,m  n 

(A16) 

16 
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E{\G  -  £(G)]2}  -  £{(G)2)  -  [£(G)]2 

=  [(!  4-  a 2)  (2)  -  <t>2(l)]  ej2n* ' 

o-,2=  'h  [1  +  cr2  —  <t» 2 ( 1 )  1  £^„2  +  '/:[(!  +  a2)  <t>  (2)  -<J>2(1)1  £a2  cos  2«m, 

n  n 

cr 2  -  Vi  [1  +  O-2  -  4>2(1)]  £/<n2  -  Vi[(l  +  tr2)  <t>  (2)  -  <D2  (1)]  ^A2  cos  2mm. 


ana 


or  I,  =  Vi  [d  +  a-2)  <t>  (2)  -  <t»2(I )]  £/l„2  sin  2 n/x. 


(A17) 

(A18) 

(A19) 


17 
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Appendix  B 

PROOF  THAT  THE  ILLUMINATION  FUNCTION  IS  BOUNDED 


Cauchy’s  inequality  is  given  by  [B 1 1 


I  v,  y, 


€ 


/-I 


Let  Y,  -  1,  then 


or 


When  T,  =  1, 


therefore 


(IT,)2  <  (I X,2)N, 


IT,2 

/ 


IT,2  _  l 
(IT,)2  A ' 


I  V,2 

(IT,)2 


is  bounded. 


REFERENCE 


Bl.  E.F.  Beckenbach  and  R.  Bellman,  Inequalities,  Springer- Verlag,  Berlin,  1971,  p.  2. 


Appendix  C 

DERIVATION  OF  VARIANCES  FOR  A  PLANAR  ARRAY 


For  the  noncorrelated  case, 

G(fi.v)  -  XX  o  +&nm)  exp  o '  exp  +  nv )]  A m,  (Cl) 

n  m 

EIGW.k)]  -  XX4*  MA*m  exp  \j(mii  +  nv)],  (C2) 

n  m 

I  Cr  (m-  .V  )  l2  —  0  +8  Ml)  (1  +8/»)  exp  L/(^/un  —  4> «)] 

n  m  r  s 

■  exp  (y[(m  -  s)/x  +  (n  -  rVl},  (C3) 

EllG^.x)!2!  “  III  Y.AnmA,&}m  0)  expi/lm  -  s)  +  (n  -  r)v]}  +  H-4i  (1  +  v2nm).  (C4) 

n  m  r  s  n  m 

n  *  r  or  m  *  s 

\E[G(fx  ,v)]  |2  —  XXXX4*™  (D^nm^rs  exp  [j[(m  -  s)n  +  (n  -  r)v]],  (C5) 

n  m  r  s 

cr  2  +  <r 2~  E[\Gin  ,v)\2\  —  \E[GifL  ,v)]\2 

-IIIl+^-^O)],  (C6) 


n  m 

[G(#i,i/)]-  XXXX  (l+8„m)  (l+8«)  exp  exp  [f(m  +  s)  n  +  (n  +  r)v],  (Cl) 

n  m  r  s 


fUG^.v)]2)  -  XXX  Xa» mArs®nm  0)  exp  [f(m  +  s)fl  +  (rt  +  r)v] 

n  m  r  s 
n  r  or  m  *  s 


t-  XX  (1  ^nm(2)A2m  exp  \J  (2mn  +  2 nv)], 

n  m 

(C8) 

[E[G(fi  ,v 

)]}2-  XXXXy4<im'4«4>J  (1)  exp  [Jim  +  s)m  +  (n  +  r)v], 

n  m  r  j 

(C9) 

<r,2  - 

a2-  XX  I(1  +<TL)*nm  (2)  -*i,0»  A„i  cos  (2 mu  +  2 nv), 

n  m 

(CIO) 

<r  12  —  Vi  XXl<1  (2)-4»i„  (1)1  aL  sin  (Imp  +  2nv), 

n  m 

(Cll) 

or  {*-  '/l  XX  IC<WI  +  Aw.  COS  (2 «M  +  2«I')] 

n  m 

(Cl  2a) 

and 

<r}“  lh  XX  iQm  -  Aw.  COS  (2 mfl  +  2nv)] 

it  m 

(Cl  2b) 

where 

c  -  (1  ^  -  <t>l„  (1)]  A2m 

(Cl  2c) 

and 

D  -  1(1  +0-2,)  (2)  -  4»2b(1)1 

19 

(Cl  2d) 

_ 

■iLf  **  V,  ^  .  ‘  *-  w  —  •-  'it 

J  K  HSIAO 


v  ■  -  • 


For  the  case  of  correlated  errors, 

G(n, v)~  Z  (1  +8«)  exP  0<*„)  Z  U  +8mB)  exp  0</>nm)  Anm  exp  L/(^M  +  mOl, 

n  m 

E[G(^,v)]~  exp  o  W  +  w)!. 

n  m 

|G(ju,i')|2  =  ZZ^  +  8„)  (1  +  8r)  exp  ~  4>r)  1  ZZ^  +  8nm)  (1  +  8„) 

nr  ms 

•  AnmArs  exp  \j(<t>nm  -  <t>r5 )]  exp  [j(m  -  s)  fi  +  («  -  /-)*/], 

£[|G(m,«')P]  -  X  Z  Z4*2*  exP  l/(™  “  SV  +  (»  ~  'M 

/if  mi 

n  ^  r  m  *  s 

+  Z  0  +  o- 2)  Z  Z4* <?»  (lM«n^«  exp  t/(m  -  s)/il  +  £  Z4*2  (1) 

n  ms  nr 

m  *  s  n  *  r 

■  J^nm^A„mArm  exp  U'(»  -  r)v\  +  ZZ  (1  +a2)  (1  +o2m)  A„2m, 

m  n  m 

|£[g<m,oii2-  ZZZZ^^O)4^  (1  >AnmA„  exp  [j[(m  -  s)n  +  (n  -  r)*]|, 

n  m  r  s 

<r?+<r i-  XIKl  +  <r*)  (1  +a„V)  -  <Pn2  (l>*i,  (1)]  AL 

n  m 

+  Z^1  +  O  -  LL^om^ns  COS  (m  -  S>M 

/!  ms 

[G (n  ,v )]2  -  ZZd.+  ®«)  (1  +8,)  exp  \j(<bn  +0,)]  ZZ*1  +  snm)  d  +  6«) 

r.  r  ms 

■  exp  \j(4>nm  +<#>„)!  AnmArs  exp  O' [Cm  +  s)n  +  (n  +  r)  x]), 
£{[G</i.»')]2}  “  Z  Z4*-  (1)  Z  Z^™  o )AnmArs  exp  {jl(m  +  s)/x  +  (n  +  rOll 

nr  ms 

n  *  r  m  *  s 

+  Z  (1  +<r2)<M2)  Z  exp  OKm  +  sV  +  2m']} 

n  ms 

m  *  s 

+  Z  Z4*2  (1)  z QlmAnmArm  exp  U  (2m#i  +  (n  +  r)v  ] 

nr  m 

n  *  r 

+  ZZ  l(1  +  0  (1  +<TnJ  <X>„(2)«l>„m(2 )]A„i  exp  L/(2 mu  +  2 nv)), 

n  m 

{£IG<m.«')]}2-  ZZZZ4*2  0)  ^L(^AnmA„  exp  { jUm  +  s)fi  +  (n  +  rV]}, 

n  v  m  j 

ZZ  l«+0  (l+o-2J«tB(2)«l>„m(2) 

n  m 

-  <f>„2(l)<I>2m(l)]  A„m  cos  (2 mu  +  Inv) 

+  ZZZ  10  +o-»)  (2)  -  ♦2mO)M>il-4mi^„ 

n  m  j 

•  cos((m  +  s)/*  +  2/iv], 


(C13) 
(Cl  4) 

(C15) 


(C16) 

(C17) 

(C18) 

(CI9) 


(C20) 

(C21) 


(C22) 
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*12-%!!  M  +0  (1  -H  ar2m)  cj>„(2)cj>„m  (2)  -  <l>2(l)«J>2m(l)l 

rt  m 

•  /ln2m  sin  (2 mu  +  2nv) 

+  Vi  HI  1(1+0*,  (2)  —  *2  (1)1 

n  m  s 

■  <t>  L  (DA^Ak  sin  l(m  +  s)  m  +  2 nv\. 

(C23) 

and 

< r,2  =  £  +  F 

(C23a) 

where 

<t\~  E  —  F. 

(C23b) 

E-'h  IL  [(1  +  0  (1 +<r2m)-<t>2(l)<I>2m(l)]  A2m 

n  m 

nm/4„  cos  (m  -  s)fi  (C23c) 

n  ms 

and 

F-  ‘AH  1(1+0  (1  +o-2m)  «I>n(2)4»„m(2)  -  Ol)*2,  (1)1  A2m  cos  (2mM  +  2np) 

n  m 

+  Vi  III  [(1  +  0  *,(2)  -  *2(l)l*i.<l)4*4«  cos  [(m  +  s)t *  +  2m>].  (C23d) 


21 
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